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Abstract 

(N ■ 

. Motivated by the recent PAMELA and ATIC results, we calculate the 

' electron and positron fluxes from the decay of lightest-superparticle (LSP) 

' dark matter. We assume that the LSP is the dominant component of dark 

I matter, and consider the case that the i?-parity is very weakly violated so that 

On , the lifetime of the LSP becomes of the order of 10^^ sec. We will see that, with 

O ■ 

. such a choice of the lifetime, the cosmic-ray electron and positron from the 

■ decay can be the source of the anomalous e fluxes observed by PAMELA and 

^ ' ATIC. We consider the possibilities that the LSP is the gravitino, the lightest 

neutralino, and scalar neutrino, and discuss how the resultant fluxes depend 
on the dark-matter model. We also discuss the fluxes of 7-ray and anti-proton, 
and show that those fluxes can be consistent with the observed value in the 
parameter region where the PAMELA and ATIC anomalies are explained. 



1 Introduction 



In astrophysics, the existence of dark matter is almost conclusive. According to the 
recent survey of WMAP pLj, it accounts for 23 % of the total energy density in the 
universe. In the standard model of particle physics, however, there does not exist 
candidate for dark matter, which is one of the reasons to call for new physics beyond 
the standard model. Supersymmetry (SUSY) is a promising model which can give 
an answer to the question; in the framework of SUSY, lightest superparticle (LSP) 
is a viable candidate for dark matter. 

The fluxes of high energy cosmic rays give information about the properties 
of dark matter. In the recent years, accuracy of the measurements of the fluxes 
have been significantly improved. In particular, anomalous signals are reported by 
PAMELA [2j and ATIC [3J in the observations of high energy cosmic-ray positron 
and electron. The PAMELA and ATIC results have attracted many attentions be- 
cause the anomalies may indicate an unconventional nature of dark matter. In fact, 
a sizable number of dark-matter models are proposed to explain the anomalies after 
the announcements of the PAMELA and ATIC results. Roughly speaking, the pos- 
sibilities to produce such high energy electron and positron are categorized into two: 



the decay and the annihilation of dark matter J^^l (For early attempts to calulate the 
spectra of cosmic-ray e±, see [3 El [71 El S M M 13 [12 [H [13 [ISl [HI 111 III EQl 
and [221 1231 [210 123 12S1 [23 IMl 123 [30] for decaying and annihilating dark matter, 
respectively.) In general, the latter case has the difficulty to reproduce the anoma- 
lous cosmic-ray positron and electron fluxes without large enhancement factor, called 



boost factor 



#2 



In the former case, on the other hand, the observed anomalies can be 
well explained with the appropriate choice of the lifetime of dark matter especially 
in leptonically decaying scenarios. 

In usual super symmetric scenario, i?-parity conservation is assumed, which pro- 
tects LSP from decaying into standard model particles and makes it a viable candi- 
date for dark matter. If we consider the case that i?-parity is violated, LSP is no 
loner stable; however, if i?-parity violation (RPV) is weak enough, the lifetime of 
the LSP can be much longer than the present age of the universe and LSP can play 
the role of dark matter [3lj. In addition, when the size of the RPV is properly cho- 



*^The other possibilities to enhance the fluxes considered using nearby pulsars [4]. 

'^^However, cosmic-ray fluxes may be enhanced without large boost factor with the Breit- 
Wigner enhancement of the annihilation cross section [31j , Sommerfelt enhancement f32' , or with 
a nearby clump of dark matter [33j . 
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sen to give the lifetime of 0(10^^ sec), produced cosmic-ray positron gives excellent 
agreement with PAMELA data |7|. 

In this paper, we calculate fluxes of cosmic-ray positron, electron, 7-ray, and 
anti-proton in various LSP dark matter scenarios, paying particular attentions to 
the results given by PAMELA and ATIC We consider the cases where the LSP 
dark matter is unstable, assuming that -R-parity is (very weakly) violated. Then, 
we compare the calculated fluxes with the results of observations. We will see that 
the PAMELA and ATIC anomalies are simultaneously explained if the lifetime of 
the LSP dark matter is 0(10^^ sec) and the mass is ~ 1 — 1.5 TeV. In addition, 
in some cases, 7-ray and anti-proton are also produced by the decay of the LSP. 
We will see that, taking account of the uncertainties in the Galaxy and propagation 
models as well as the error in the observations, the scenarios are not excluded by the 
observations of 7-ray and anti-proton fluxes yet. 

The organization of this paper is as follows. In the next Section, we explain the 
basic procedures to calculate the high energy cosmic-ray fluxes. Then, we discuss the 
cosmic-ray fluxes for the cases where the LSP is the gravitino, the lightest neutralino, 
and the sneutrino in Sections [31 HI and 0, respectively. Section [H] is devoted for 
conclusions and discussion. 

2 Basic Setup 

As we have discussed in introduction, we consider the case where the LSP is unstable 
and has lifetime much longer than the present age of the universe so that most of the 
LSPs produced in the early universe survive until today. Then, if the relic density 
of the LSP is right amount, the LSP can be dark matter. In such a case, the LSP 
dark matter becomes the source of high energy cosmic rays. 

The fluxes of the cosmic rays from the decay of the relic LSP depends on what 
the LSP is, and how the LSP decays. The LSP should be an electrically neutral 
particle and we consider the following three important cases: 

• Gravitino LSP 

• Neutralino LSP 

• Sneutrino LSP 

In addition, we adopt -R-parity violation so that the LSP becomes unstable. Then, 
the lifetime of the LSP becomes longer as the coupling constants for the RPV interac- 
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tions become smaller. We therefore treat the lifetime of the LSP as a free parameter 
in the following analysis. To be specific, we assume that the i?-parity violation is 
so weak that the lifetime of the LSP becomes much longer than the present cosmic 
time. 



2.1 Sources of the Cosmic Rays 

With the decay of relic LSP, energetic particles (in particular, 7, e^, and anti-proton 
p) are produced. The production rate of the energetic particle / is given by 



n.sp 



dNi 

He 



(2.1) 

- dec 

where [dN i / dE]^^^ is the energy distribution of particle I from the single LSP decay 
process. We have used the PYTHIA package |35j to calculate [dNi / dE]^^^. In 
addition, tlsp and nLSP are the lifetime and number density of the LSP at the 
present universe, respectively. In the calculation of the cosmic-ray spectrum, the 
origin of the dark-matter LSP is unimportant. Even though it is often assumed that 
the relic density of the LSP is thermally determined, non-thermal production of the 
LSP dark matter is also possible [SB]- We thus do not specify the origin of the relic 
LSP in the following analysis, and set nLSP to be I^lsp = ^dm, where II^sp arid VL^m 
are the density parameters of the LSP and dark matter, respectively!^ Then, the 
mass density of dark matter is given by pDM(a^) = '^LSP'T-Lsp(a^)- In calculating the 
fluxes of high-energy cosmic rays, we adopt the Navarro-Frank- White (NFW) mass 
density profile pTTj : 

PNFw(x)=Pe '^^t'^r , (2.2) 

where Pq ~ 0.30 GeV/cm^ is the local halo density around the solar system, Tc — 
20 kpc is the core radius of the dark matter profile, Tq ~ 8.5 kpc is the distance 
between the Galactic center and the solar system, and r is the distance from the 
Galactic center. Using the Qj given in Eq. (12. ip as a source term, we solve the 
propagation equations for individual particles. The diffusion zone is approximated 
as a cylinder with half-height L and radius i? = 20 kpc. 



^'^Even if Olsp < ^^dm, the relic LSP can still be the source of high energy cosmic rays. Then, 
the fluxes can be obtained by rescaling the lifetime as tlsp — > (f^LSp/^DM)''"LSP- 
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Ml 


MED 


M2 


L [kpc] 


15 


4 


1 


5 


0.46 


0.70 


0.55 



Kf2 [kpcVMyr] 0.0765 0.0112 0.00595 
Table 1: Parameter sets for the propagation model of e 



2.2 Electron &; Positron fluxes 

In order to calculate the fluxes of electron and positron from the LSP decay, we 
derive a static solution of the following diffusion equation: 

= K,^{E)V'U{E,x) + ^ [h{E)U{E,x)]+Q,^{E,x), (2.3) 

with the condition /e± = at the boundary of the diffusion zone, where /e± is the 
number density of per unit energy. Our basic procedure to solve the diffusion 
equation (12.31) is explained in Appendix |Al 
We approximate the function Kf.± as [38j 

= irfi^SeV, (2.4) 
with -E'ceV being energy in units of GeV, while the energy-loss rate h is given by 

6 = 1.0 X lO^i^E^^v GeV/sec. (2.5) 

We adopt three sets of diffusion parameters, which are summarized in Table [U The 
MED set gives the best-flt value in the boron-to-carbon ratio (B/C) analysis, while 
the maximal and minimal positron fractions are expected in the Ml and M2 sets 
without conflicting the B/C analysis. 

Once /e± are given by solving the above equation, the fluxes can be obtained as 

where Xq is the location of the solar system, and c is the speed of light. In order to 
calculate the total fluxes of e^, we also have to estimate the background fluxes. In our 
study, we adopt the following fluxes for cosmic-ray electrons and positrons produced 
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by collisions between primary protons and interstellar medium in our galaxy |39j : 



-1.1 



[$e-]™r. = ' o 9TT(GeV cm^ sec str)-\ (2.7) 

["^e-Lr = ^ T^(GeV cm^ sec str)-\ (2.8) 

[$e+],.. = ^'^^ ^(GeV cm2 sec str)-^ (2.9) 

^ ' ^'^^^ 1 + 650^2-3^ + ISOOE^fv ^ ^ ^ ^ 

With these backgrounds, the total fluxes are obtained as 

i^e-U = [^e-]^M + [^e-U^+[^e-L,. (2.11) 

Using the fluxes deflned above, the positron fraction, which is measured by the 
PAMELA, is deflned as 

p ^ [$e+(^)]tot .2 12^ 

[<l>,-(E)]tot + [<fe+(i?)]tot' ^ ■ ' 

2.3 Anti-proton flux 

The flux of p from the LSP decay is obtained by solving the diffusion equation: 

^MEl^ = K,{E)V'U{E,x)-^^[V3\gn{z)U{E,x)] 

-2h6iz)T^^Jp{E,x) + QpiE.x), (2.13) 

where z is the distance from the Galactic plane. Here, Vc is the convection velocity, 
h is the half height of the thin Galactic disc, which is taken to he h = 100 pc, and 
Tann IS the annihilation rate of p in the Galactic disc, which is given by 

Tann = (^H + 4^^^niie)(TppVp, (2.14) 

where we use the number density of Hydrogen and Helium in the Galactic disc to 
be nn = 1 cm"^ and nue = O.OTnn- The cross section app is given by [iOl 1^ 

. 661(1 + O.OllST^eV - 0.948TO-o^5i) mb : Tcev < 14.6 GeV ^ , 
' 36T-°v^ mb : Tcev > 14.6 GeV 

with Tcev being the kinetic energy of the anti-proton in units of GeV. In addition, 
as in the case of e^, the function Kp is parametrized as 

Kp = K^'^PppLv, (2.16) 



MAX MED MIN 



L [kpc] 15 4 1 

5 0.46 0.70 0.85 

Kf^ [kpcVMyr] 0.0765 0.0112 0.0016 

K [km/s] 5 12 13.5 

Table 2: Parameter sets for the propagation models of p. 



where (3p is the velocity of the anti-proton, and pcev is the momentum in units of 
GeV. Parameter sets for the diffusion equation used in our analysis are summarized 
in Table [2l Again, the MED set gives the best-fit to the B/C analysis, while the 
maximal and minimal anti-proton fluxes are expected in the MAX and MIN sets. 
Once fp is obtained, the anti-proton flux at the solar system is calculated as 

[%iE%M='%fv{E,S^). (2.17) 



2.4 7-ray flux 

The flux of the 7-ray is calculated by the sum of two contributions: 

['^'7]dM = ['^'7]™ + ['^'Tlhalo , (2-18) 

where the flrst and second terms in the right-hand side are fluxes of 7-ray from 



cosmological distance and that from the Milky Way halo, respectively 
The flux from cosmological distance [$^1 is obtained as 

<^ L yJcosmo 



#4 



1 /""^ 

[%lo.mo = -—-— dE'G,{E,E') 

"^LSP^LSP Je 



Ie' 



(2.19) 



dec 



Here, the propagation function of 7-ray turns out to be 



*'*7-ray may be also produced by the inverse Compton (IC) scattering process. We have estimated 
the 7-ray flux from the IC process with the cosmic microwave background radiation in the sky 
region 10° < b < 20° (with h here being Galactic longitude) to compare with the preliminary 
FERMI results, and found that the flux is much smaller than the FERMI data. Detailed study of 
the 7-ray flux from the IC process will be given elsewhere [43] . 
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where Hq is present Hubble expansion rate, pc is critical density, while Q.m 
0.137 /i"^ and Q/^ ~ 0.721 (with h__^ 0.701) are density parameters of total matter 
and dark energy, respectively [1]. 
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On the other hand, the flux from the Milky Way halo [^'7]haio calculated as 



1 

["^7] halo ~ 



I PDM{l)dl) , (2.21) 

dec \>^l.o.s. / dir 



where the integration should be understood to extend over the line of sight (l.o.s.) 
and (■■■)dir means averaging over the direction. In the EGRET observation [12] . 
the signal from the Galactic disc is excluded in order to avoid the noise. Thus, in 
order to compare our numerical results with the EGRET observation, we exclude 
the region within ±10° around the Galactic disc in the averaging. 

For the background flux against the signal, we adopt the following flux formula 
which is estimated from the EGRET observation in the energy range 0.05 GeV < 
E < 0.15 GeV mi: 

E^ [^'tIbg - 5-18 X 10"^ (cm^ sec str)"^ GeV x ( -^-^ j . (2.22) 

We have assumed that the spectrum of the background flux from astrophysical origins 
follows a power law, and its behavior can be extracted to the high energy region. 
The total 7-ray spectrum is then given by 

+ [%]bg- (2-23) 



3 Gravitino LSP 

Now, we are at the position to discuss the fluxes of for individual dark matter 
scenarios. The first example is the case where the gravitino, which is denoted as ip^j,, 
is the LSP and hence is dark matter. The model discussed here is the same as that 
given in [331 |35l |36] , and we consider the following RPV interaction 

Cjipv = BiLiHu + m\^^liHl + h.c, (3.1) 

where Li is left-handed slepton doublet in i-th generation, while Hu and Hd are up- 
and down-type Higgs boson doublets, respectively. (In the present analysis, H^, and 

We found a typo in Eq.(4.3) of [Hj. in the first factor of right-hand side (i.e., 

1 /2 

'^T^Hf)Q.M) should be fij^j . 
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Li are defined in the frame in wliicli tlie bi-linear i?-parity violating terms in tfie 
superpotential vanisli.) 

Tlie effects of tlie RPV with (13.11) is parametrized by the VEV of the sneutrino 
fields z/j. To parametrize the VEV of z/j, we define 

ly) Bisin 13 + m] cos/? 
Ki = ^ = , (3.2) 

V m~ 

where v ~ 174 GeV is the VEV of standard-model-like Higgs boson, tan/5 = 
{H^)/{H^), and mp- is the mass of Since we are interested in the case of very 
long lifetime, we consider the case that <^ 1. 

With the RPV operator given in Eq. (13.11) . the gravitino decays as -0^ — *• jUi, 
Zui, Wli, and hui. (For detailed calculations of the decay rates for these processes, 
see in [S].) In particular, in the limit that the gravitino is much heavier than the 
weak bosons, the following relation holds: T^^^zu, — ^■4>f,^hui — ^rv-p^iy/,! and 
the process ip^^ — > 7z/j is suppressed. Then, once the relic gravitino decays, the 
produced /j (as well as the weak and Higgs bosons) becomes the source of cosmic-ray 
electron and positron. In addition, with the hadronic decay of the weak and Higgs 
bosons, energetic anti-proton and 7-ray are also produced. Fluxes of these particles 
have been measured, and in the following, we compare the expected fluxes with the 
results of observations. 

As we have seen, the fluxes of the cosmic rays originating from the gravitino 
decay is inversely proportional to the lifetime of the gravitino r3/2. With the RPV 
interaction given in Eq. (13. ip . the lifetime of the gravitino is approximately given by 

.3,.«xl0^=.cx(^)-(^)- ,3.3, 

where 

(3-4) 



Now, we show the cosmic-ray fluxes originating from the gravitino decay. First, 
we consider the fluxes of e^, motivated by their observed anomalous fluxes recently 
reported by the PAMELA and ATIC experiments. In order to discuss the preferred 
lifetime to explain the anomalous positron fraction observed by the PAMELA ex- 
periment, we define the variable as 

/p(obs) _ p(th)\2 

^ =2^ jM-2 ' (3.5) 



where R °+ ^ and ■ are positron fraction in i-th bin measured by the PAMELA 
and that predicted in the unstable dark matter scenario, respectively, and SR^^^f is 
the error in the observed fraction. Since the positron flux in the low-energy region 
is sensitive to the background fluxes, we only use the data points with _E > 15 GeV 
(5 data points) in the calculation of x^- (As we will discuss, the positron fraction 
depends on the background. We use the variable just to estimate the preferred 
value of the lifetime to explain the PAMELA results.) 

The best-fit lifetime to explain the PAMELA anomaly depends on the gravitino 
mass, the flavor of the final-state lepton, and the propagation parameters. In our 
analysis, for simplicity, we consider the cases where the LSP dominantly decays into 
fermions in one of the three generations. Concerning the other parameters, we will 
discuss how the resultant fluxes depend on them. 

For the case where the gravitino decays only into the first-generation lepton, we 
plot the positron fraction in Fig. [T]with the best-fit lifetime. Here, we show results 
with adopting the MED and M2 propagation models, because the results with the 
Ml and MED models are almost the same. For the MED (M2) propagation model, 
the best-fit lifetime is given by 2.0 x 10^^ sec, 1.1 x 10^^ sec, and 8.6 x 10^^ sec 
(9.3 X 10^5 sec, 5.0 x 10^^ sec, and 4.3 x lO^^ sec) for mg/a = 300 GeV, 600 GeV, 
and 1.2 TeV, respectively. 

The positron fraction for the cases where the gravitino decays only into second- 
and third-generation lepton are shown in Figs. [2] and [3l respectively. Here, we show 
the results with the propagation models MED and M2. We can see that the predicted 
positron fraction well agrees with the PAMELA result if the lifetime is properly cho- 
sen. In particular, when the final-state lepton is first- (second-) generation, the fit 
is excellent irrespective of the gravitino mass with MED (M2) propagation model. 
We note here that, for e~, the background flux is significantly larger than the signal 
flux, while the dark matter contribution dominates for e"*". Thus, the result is sen- 
sitive to the choice of background because we plot the positron fraction. However, 
the theoretical calculation of the positron fraction contains parameters both in the 
particle-physics model (i.e., the lifetime and the flavor of the final-state leptons) and 
in the propagation model. Thus, we believe that positron fraction observed by the 
PAMELA can be explained in the present scenario with other choice of the back- 
ground fluxes. For example, even if the normalization of the background e~ flux is 
changed, we can obtain almost the same positron fraction by varying the lifetime. 

Next, we consider the total flux $e+ + "^e-- The numerical results are shown 
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Figure 1: Positron fractions in (a) MED and (b) M2 models for the case where the 
gravitino dominantly decays to the first-generation lepton. Dot-dashed Une is the positron 
fraction calculated only by the background fluxes. Here, we take ms/2 = 300 GeV, 600 
GeV, and 1.2 TeV (from left to right) with 2.0 x 10^^ sec, 1.1 x 10^^ sec, and 8.6 x 10^^ sec 
(9.3 X 10^^ sec, 5.0 x 10^^ sec, and 4.3 x 10^^ sec) in MED (M2) model, respectively. 
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Figure 2: Same as Fig. [U but for the case where the gravitino dominantly decays to 
second-generation lepton. Here, we take r3/2 = 1-5 x 10^^ sec, 1.1 x 10^^ sec, and 8.6 x 10^^ 
sec (9.3 X 10^6 sec, 5.8 x lO^^ sec, and 5.0 x lO^^ sec) in MED (M2) model, which are 
the best-fit lifetime. 
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Figure 3: Same as Fig. \T5\ but for the case where the gravitino dominantly decays to 
third-generation lepton. Here, we take = 9.3 x 10^^ sec, 8.6 x 10^^ sec, and 7.9 x 10^^ 
sec (6.3 X 10^^ sec, 5.4 x 10^^ sec, and 5.4 x 10^^ sec) in MED (M2) model, which are 
the best-fit hfetime. 
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Figure 4: Total fluxes of positron and electron in (a) MED and (b) M2 models for the 
case where the gravitino dominantly decays to the first-generation lepton. Dot-dashed line 
is the background flux. Here, we take the the same mass and lifetime as Fig. [H and also 
plot PPB-BETS data g7]. 
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Figure 5: Same as Fig. |H but the gravitino dominantly decays to second-generation 
lepton. We take the the same mass and hfetime as Fig. [21 and also take 01,3/2 = 2 TeV 
with = 7.4 X 10^^ sec and 4.6 x 10^^ sec in (a) and (b), respectively. 
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Data: ATIC,PPB-BETS 
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Figure 6: Same as Fig. HI but the gravitino dominantly decays to third-generation lepton. 
We take the the same mass and hfetime as Fig. [3l 
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in Figs, m — [6] for the cases where the gravitino decays only into first-, second-, or 
third-generation lepton. For the calculation of $e+ + "^e- ? we use the best-fit lifetime 
to fit the PAMELA data. In the figures, it can be seen that the observed anomalous 
structure is well reproduced with both MED and M2 models by relevantly choosing 
m3/2, except for the case that final-state lepton is third-generation. For the case 
of final-state lepton being in the first- (second-) generation, the result is a good 
agreement with the observation when m3/2 ~ 1.2 TeV (2 TeV). Here, we note that 
the total flux is not sensitive to the background because the signal from the dark 
matter is larger than (or at least comparable to) the background. 

With the decay of the gravitino dark matter, energetic 7-ray and anti-proton are 
also produced. Thus, with the observations of the fluxes of these particles, we may 



confirm or exclude the present scenario Notice that 7 and p are mostly from the 



hadronic decays of the weak and Higgs bosons. Thus, the fluxes of these particles are 
insensitive to the flavor of the final-state lepton, except for the case where r-lepton 
is produced by the decay. (See the following discussion.) 

In Fig. [3, we plot the flux of the 7-ray in the present scenario with the primary 
lepton being in flrst- or second-generation. From the flgure, we see that the expected 
7-ray flux well agrees with EGRET data irrespective of Because the flux from 

the gravitino decay is larger than the background, the 7-ray flux does not highly 



depend on the background]^ On the other hand, if the primary lepton produced in 
the decay is in third-genaration, 7 from ir^ decay also contributes to the total flux 
jSlj . In Fig. [S], we show the result for such a case. From the flgure, one can see an 
increase of the flux in the high energy region of E > 100 GeV. 

The anti-proton flux from the decay of the gravitino dark matter is shown in 
Figs. M — [m We can see that the anti-proton flux depends on the propagation 
model. However, with the MED and MIN models, the flux from the gravitino decay 



#^High energy 7-ray from the Galactic center has been calculated in the scenario where dark 
matter annihilates into W^W^ pair [151 [TH]. We have checked that the high energy 7-ray flux 
in the present scenario is much smaller than that in the annihilating scenario (into W~^W~ pair), 
assuming that the PAMELA and ATIC anomalies are from the decay or the annihilation of dark 
matter. Then, we found that the total 7-ray flux in the present scenario is consistent with the 
HESS observation g^. 

'^^The recent preliminary results of the FERMI satellite (for the region 10° < b < 20°) indicates 
no anomalous excess in the high energy 7-ray ^.50] • We have also calculated the 7-ray flux for the 
region 10° < b < 20° in the present scenario, and found that the dark-matter contribution to the 
flux for such a region is smaller than the FERMI data. 
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Figure 7: 7-ray flux. Dot-dashed line is the background flux. Here, we take mass and 
lifetime as the same as figure (a) in Fig. [T] 
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Thus, we conclude that 



is comparable to or smaller than the observed fluxes, 
the present scenario is not excluded by the observation of the cosmic-ray p flux, 
taking account of the uncertainties in the propagation model and estimation of the 
background. However, if a better understanding of the propagation of the anti- 
proton becomes available in the future, it will provide a significant test of the present 
scenario. 

Before closing this section, we comment that this scenario may be tested by the 
LHC experiment. Indeed, in this scenario, the lightest superparticle in the MSSM 
sector (which we call MSSM-LSP) may decay inside the detector. In the present 
scenario, the MSSM-LSP is likely to decay to the standard model particles via the 
RPV interaction even though there exists a superparticle (i.e., gravitino) lighter than 
the MSSM-LSP. As we have mentioned, the n parameter is expected to be O(10~^°), 
which gives the lifetime of the MSSM-LSP of the order of ~ io~(^~^) sec. Thus, 
the typical decay length of the MSSM-LSP is expected to be much longer than the 
sizes of the ATLAS and CMS detectors. However, since enormous amount of SUSY 



We also estimated anti-proton to proton ratio to compare recent data by PAMELA [55] , using 
the proton flux observed by BESS [53] and CAPRICE [54|. Then we found that the p/p ratio is 
order of magnitude smaller than the PAMELA data if we take the MIN propagation model and that 
it is comparable to or a few times larger than the PAMELA data with MED model. Thus, taking 
account account of the uncertainties in the propagation model as well as those in the background 
proton flux, the present scenario is not excluded yet. 
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Figure 8: The same as Fig. [71 but for the case where primary lepton produced in the 
decay is third-genaration. 



events is expected at the LHC experiment, some of the MSSM-LSP produced at the 
LHC may decay inside the detector, which resuhs in drastic signal, hke a displaced 
vertex or a kink in a high-energy charged track. In addition, if such decay processes 
can be observed, it may be also possible to constrain the lifetime of the MSSM-LSP, 
which may be used for the determination of the k parameter [57] . 



4 Neutralino LSP 

In the previous section, we have considered the case that the gravitino is the LSP. In 
conventional SUSY models, another important candidate for the LSP is the lightest 
neutralino, which is a linear combination of Bino, neutral Wino, and two neutral 
Higgsinos. To make our discussion simple, in this Section, we assume that the lightest 
neutralino is (almost) Bino-like. This is the case, for example, in large fraction of 
the parameter space of models with the grand-unification condition among gaugino 
masses. 

With the RPV operators given in Eq. fl3.1l) . the Bino-like neutralino B decays as 
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Figure 9: Anti-proton flux in MIN, MED, and MAX models. Here, we take 01,3/2 
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2.0 X 10^^ sec, which is the best-fit lifetime in the case that the 



gravitino dominantly decays to first-generation lepton, and also plot the observation data 
by BESS [55] and CAPRICE [56]. 



B Zui, Wli, and hui. Decay rates for these processes are given byc^ 
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(4.1) 
(4.2) 
(4.3) 



B , 



where gz = ^/ gj + gi (with gi and g2 being the gauge coupling constants of the 
f/(l)y and SU{2)l gauge groups, respectively), 6y/ is the Weinberg angle, is 
Bino-like neutralino mass, and m^, miy, Tnh are masses of corresponding gauge and 
Higgs bosons. 

As one can see from the above decay rates, we obtain the relation T j^^^u — 
^B^hi^i — ^^B-^wk- Remember that, for the case of the gravitino LSP, the same 
(approximated) relation holds. Thus, the fluxes of the high energy cosmic rays in 
the Bino LSP case is expected to be similar to that in the gravitino LSP case as 
far as the lifetime is ~ 10^^ sec. Since the decay rate of the Bino is not suppressed 



#^Here, we assume that the lightest Higgs boson is almost standard-model like. 
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Figure 10: Same as Fig. [9l except for taking 7713/2 = 600 GeV and corresponding best-fit 
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Figure 11: Same as Fig. [9l except for taking 7773/2 = 1-2 TeV and corresponding best-fit 
lifetime. 
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by the Planck scale, the size of the ki parameter relevant to explain the PAMELA 
and ATIC anomalies is much smaller than that in the gravitino LSP case. Indeed, 
with the decay rates given in Eq. (14.11) — (14.31) . the lifetime of the Bino LSP in the 
present case is estimated as 



r^^2xl025secx (^rV^^) . (4.4) 

So far, we have considered the bi-linear RPV interaction given in Eq. (13.11) . In 
such a case, as we have discussed, high energy 7 and anti-proton are also produced 
by the decay of the LSP. However, if we consider other types of RPV interactions, 
it may be possible to enhance the e"*" and e~ fluxes without affecting 7-ray and 
anti-proton fluxes. This is the case where the Bino LSP decays mainly via the PRV 
superpotential: 

W^RPV = IXijkULjEl, (4.5) 

where Xijk = —Xjik- With this superpotential, the Bino decays as B —> i^dL,j^R,k 
and i^jl'iit^k diagrams with a (virtual) slepton propagation. Hereafter, let us 
consider the e"*" and e~ fluxes in such a case. 

With the above superpotential, the Bino decays into the three-body final state, 
and hence the final-state leptons are not monochromatic. For simplicity, we consider 
the case that the right-handed sleptons are lighter than left-handed ones, so that the 
diagram with the propagator of the right-handed slepton dominantly contributes to 
the Bino decay. Then, denoting the energies of if^- and /^^ in the rest frame of B 
as Eij^ and Ei^, respectively, the energy distribution of the charged leptons for the 
process B — >• Vilf^ jVfi ^ given by 



dEiJEi^ 647r3m^ [(mr Jm^Y - 1 + Zi^_ 



2 ■ 



(4.6) 



Where 2;/^ ^ = 2Ei^ rI"^b- (Notice that < ^ < 1.) 

In Fig. [T2I we plot the energy distributions of the final-state leptons: 



dzL,R ~ T B^^^i± IT J '""'^ dzijzi^ 



L,j R.k 

Notice that this quantity depends only on the ratio mi^^/nij^. (In the figure, we 
take mj^/m^ = 1.2.) As one can see, the left-handed lepton emitted in the decay is 
likely to be energetic. Thus, if the coupling constant Xijk is sizable so that the Bino 
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3 



Neutralino LSP 




Figure 12: Distribution of the left- (solid) and right-handed (dashed) charged leptons 
emitted from the Bino decay via the LLE'^-type superpotential. Here, we take mj^^/mj^ = 
1.2. 

decays dominantly as 5 — >• e^z/j/^^, significant effects on the electron and positron 
fluxes are expected. 

In Fig. [131 we show the positron fraction in the Bino dark matter case with the 
RPV superpotential given in Eq. (14.51) . Here, we consider the case that A122 is the 
largest so that the Bino decays via this coupling, and we use mj^/m^ = 1.2. The 
best-fit lifetime to explain the PAMELA anomaly depends on as well as on the 
propagation model as in the case of the gravitino LSP. For the MED (M2) model of 
propagation, the best-fit lifetime is given by T3/2 = 3.2 x 10^^ sec, 2.0 x 10^^ sec, and 
1.2 X 10^5 sec (2.0 x lO^^ sec, 1.1 x 10^^ sec, and 5.0 x 10^^ sec) for = 300 GeV, 
600 GeV, and 1.2 TeV, respectively. In addition, in Fig.lHl we plot the flux of e~^+e~. 
We can see that the simultaneous explanation of the PAMELA and ATIC/PPB- 
BETS anomalies may be possible in the present scenario if ^ 1.5 TeV and 

~ 9.3 X lO^^sec (4.0 x lO^^ sec) for MED (M2) model. 

So far, we have considered the case that A122 is the largest. We have checked 
that similar results for the positron fraction and the total flux are obtained as far 
as the first-generation lepton is emitted in the decay. If B decays only into second- 
and third-generation leptons, on the contrary, the total flux $e+ + $e- is suppressed 
compared to the observed values when we adopt the best-fit lifetime to explain the 
PAMELA anomaly. 
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Figure 13: Positron fractions in (a) MED and (b) M2 models for the case wliere only A122 
is non-zero. Here, we take ruf^ = 300 GeV, 600 GeV, and 1.2 TeV from left to right and 
= 3.2 X 10^6 sec, 2.0 x 10^^ sec, and 1.2 x 10^6 sec (2.0 x 10^^ sec, 1.1 x 10^6 sec, and 
5.0 X 10^5 sec) in MED (M2) model, respectively. 
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Figure 14: Total fluxes in (a) MED and (b) M2 models for the case where only A122 is 
non-zero. We take the same values for m p, and Tp, in Fig. [T3l and also take 
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with TB = 9.3 X 10^5 sec (4.0 x 10^^ sec) in (a)((b)). 
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5 Sneutrino LSP 



The third candidate for the LSP is the sneutrino. In the framework of the MSSM, 
the sneutrino may not be a popular candidate for the LSP. However, the hghtest 
sneutrino can be dark matter without conflicting phenomenological constraints, as 
we discuss below. In the previous section, we have seen that, with the LLE'^ type 
RPV superpotential given in Eq. (14.51) . enhancement of electron and positron fluxes 
is possible without affecting the 7-ray and anti-proton fluxes. If the Bino is the LSP, 
the final-state leptons are not monochromatic, so that the end-point behavior of the 
cosmic-ray is smoothed. On the other hand, if the sneutrino is the LSP and hence 
is dark matter, and also if it decays via the LLE^ type RPV superpotential, shape 
of the spectrum drastically changes and a sharp edge at the end-point can be 
obtained [7j. 

In the MSSM, there only exist three left-handed sneutrinos v^s ("L" means left- 
handed in this Section), and one of them may be the lightest superparticle. By 
tuning the soft SUSY-breaking scalar and gaugino masses, the sneutrino becomes 
the LSP in some parameter region. In addition, the direct detection constraint on 
the sneutrino dark matter [5B] can be avoided by introducing small lepton-number 
violating operator ~ {LHuY [59j. In addition, there is another possibility to realize 
the sneutrino dark matter scenario, in which the right-handed sneutrino vr becomes 
the LSP [nn]- Since various neutrino-oscillation experiments suggest that the neu- 
trinos are massive, we expect the existence of the right-handed (s) neutrinos. If the 
neutrino masses are Dirac-type, vn becomes as light as other MSSM superparticles 
in the gravity-mediation SUSY breaking scenario. (In other cases, vr may be much 
lighter than MSSM superparticles.) Thus, in some parameter space, can be the 
LSP. In such a case, can also be dark matter. The following arguments holds if 
vl or is the LSP. 

If the sneutrino in i-th generation is the LSP, it decays as z/j l^l^ assuming 
that the dominant RPV interaction is given by the LLE'^ type RPV superpotential 
given in Eq. (14.51) . The decay rate of this process is 

where 6^, is the mixing parameter in the sneutrino sector. When the vl is the LSP, 
Oi, = 1. Even if z/j is right-handed, it decays via the RPV interaction given in Eq. 
(14.51) because there should exist left-right mixing term of the sneutrinos, which is of 
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the form 



Clr = 6ml^-^ULi^ji + h.c. (5.2) 
In such a case, 9c, in Eq. (15. ip should be 

e, = 4^^^. (5.3) 

mf, — m- 

Thus, in any case, the sneutrino LSP decays into charged-lepton pair via the RPV 
interaction given in Eq. fl4.5p . Thus, if the sneutrino is dark matter, those leptons 
become the source of cosmic-ray e^. Importantly, in such a case, dark matter decays 
only leptonically, and the productions of hadrons and photon are irrelevant. Thus, 
the present model may produce significant amount of in cosmic ray without 
affecting the anti-proton and 7-ray fluxes. 

In the following, we calculate the electron and positron fluxes from the decay of 
sneutrino dark matter. Here, we take mc, = 300 GeV, 600 GeV, and 1.2 TeV. The 
resultant flux depends on the flavors of the final-state leptons. Here, for simplicity, 
we consider three simple decay modes: z/ e+e", fJ'~^fi~, and r+r". 

The results of the calculations of positron fraction are shown in Fig. [T5] — fTTl 
where MED mode of the propagation is adopted. In addition, we use the best-fit 
lifetime for each mp. We see good agreements with the observation for all the cases. 
We have also checked that a reasonable agreement between the prediction and the 
PAMELA data is obtained even with the M2 model if the sneutrino decays into 
fi^fi~ and T~^T~ . 

Next, we discuss the total flux $e+ + '^'e-- The numerical results are shown in 
Figs. [18] — [20I taking the same parameters as in Fig. [15] — [TJl respectively. One can 
see the anomalous behavior is explained when ~ 1.2 — 1.5 TeV in the cases of 
final state: e~^e~ and fj^'^fi~ with MED model. In addition, we checked that the same 
anomalous behavior is obtained in the case of final state: fi~^fi~ with M2 model. 



6 Conclusions and Discussion 

In this paper, we have studied the cosmic-ray fluxes from the decay of LSP dark 
matter, motivated by the recently reported anomalies by PAMELA and ATIC We 
have introduced several types RPV operators so that the LSP becomes unstable, 
and calculated the fluxes of e^, as well as those of p, and 7-ray, assuming that the 
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Figure 15: Positron fractions in MED model for the case where the sneutrino decays to 
the final state: e^e~. Here, we take = 300 GeV, 600 GeV, and 1.2 TeV from left to 
right and = 5.4 x 10^^ sec, 2.5 x 10^^ sec, and 1.6 x 10^^ sec. 
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Figure 16: Same as Fig. [15] in the case of final state: /i^/x . Here, we take Ty = 3.7 x 
10^^ sec, 2.3 X 10^^ sec, and 1.2 x 10^^ sec, which are the best-fit hfetime with MED model. 
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Figure 17: Same as Fig. [15] in the case of final state: r+r . Here, we take Ti, = 1.8 x 

1026 

sec, 1.5 X 10^^ sec, and 1.1 x 10^^ sec, which are the best-fit hfetime with MED model. 



LSP is the dominant component of dark matter. The detailed shape of the spectra 
of cosmic-ray depend on the properties of the LSP dark matter. However, when 
the lifetime of the LSP is of 0(10^^ sec), the predicted positron fraction can be 
consistent with the observed one by PAMELA irrespective of the mass of the LSP. 
On the contrary, in order to explain the ATIC anomaly for the total flux $e+ + '^'e- ? 
the mass of the LSP is required to be 1 — 1.5 TeV, depending on the decay modes of 
the LSP. If a monochromatic is produced by the two-body decay process, the LSP 
mass can be as light as ~ 1 TeV, while larger mass is required for the ATIC anomaly 
for other cases. In any case, in order to explain the ATIC anomaly, relatively large 
value of the MSSM particle masses are required, which may make it difficult to solve 
the naturalness problem with supersymmetry. 

When the LSP decays via the bi-linear RPV interaction given in Eq. (13. ip . sig- 
niflcant amount of 7 and anti-proton are also produced by the decay. In particular, 
it has been discussed that the anti-proton flux may give a stringent constraint on 
the decaying dark matter scenario [16]. However, taking account of the uncertain- 
ties in the propagation models as well as the errors in the observed fluxes, we have 
shown that the scenarios discussed in this paper are not excluded by the present ob- 
servations. In addition, we have also seen that the 7-ray flux is consistent with the 
currently available observational results. With improved knowledges about the prop- 
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Figure 18: Total fluxes of positron and electron in MED model for the case sneutrino 
decays to the final state: e~^e~ . Here, we take the the same mass and lifetime as Fig. [TC} 
and also take ma = 1.5 TeV with 1.2 x 10^^ sec. 

agation of the anti-proton, more detailed test of the scenario may become possible 
in the future. It is also notable that a precise measurement of the cosmic-ray 7 flux 
is expected by the Fermi Telescope. Furthermore, informations about the decaying 
dark matter may be imprinted in the synchrotron radiation from the Galactic center 
and in high-energy neutrino flux [HS]- Future improvements of the observations 
of cosmic-ray fluxes should provide better understandings of the properties of dark 
matter. 
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A Green's Function 

We discuss here how we obtain the Green's function for cosmic ray positrons from 
the LSP decay. As mentioned in Section [21 the Green's function is obtained by 
solving the diffusion equation (12.31) with the boundary condition fe±{E,x) = at 
the surface of the diffusion zone. Because of the condition, it is convenient to expand 
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Figure 19: Same as Fig. [18] for the case of final state: ^"'"/z . We take the the same mass 
and lifetime as Fig. [T6l 
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Figure 20: Same as Fig. [18] for the case of final state: r+r . We take the the same mass 
and lifetime as Fig. [17] 
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the solution by Bessel series for r (the radius of the cyhnder) and Fourier series for 
the z (the thickness of the zone). We then find that the positron flux is obtained as 



- - g 



dE'G{E, E') 



~dE' 



(A.l) 
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Defining the the variable 
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(A.2) 



the Green's function G{E, E') turns out to be 
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where r = E'^/b{E), Jn is the n-th order Bessel function, and Cn are successive zeros 
of Jq. The above formula is useful to calculate the Green's function, except for the 
region E ~ E'. This is due to the exponential suppression in Eq. flA.Sp . For the 
other energy region, we have checked that the summations ove r n and m are well 
converged when we take the Bessel (Fourier) series large enough 
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When E c::^ E', the convergence becomes worse. In fact, the result is not converged 
for the cuspy profiles such as NFW and Moore profiles even if we take the summations 
up to (9(1000). However, note that the flux for ~ £" originates from the LSP 
decay in the vicinity of the solar system. Using the fact, we can express the Green's 
function G{E, E') in another form as 



G{E,E') = [exp (-ir2^rX(E,E')V') pdm( 



(A.5) 



r=rQ 



In contrast with the previous formula, this formula is valid in the region X{E, E') ^ 
1, namely E c:^ E'. By matching two formulas in Eqs. (1A.3P and (lA.Sp at appropriate 
value of X{E,E'), we can obtain the Green's function numerically without wasting 
time for computation. 

#iOWe have taken the summations up to n,m — 100 in our study. 
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